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Abstract: A recently proposed form of dual theory for three dimensional superconductor 
is rederived starting from the lattice electrodynamics and studied by renormalization group. 
The superfiuid density below and close to the transition vanishes as inverse of the correlation 
length for the disorder field. The corresponding universal amplitude is given by the fixed 
point value of the dual charge, and it is calculated to the leading order. The continuum 
dual theory predicts the divergence of the magnetic field penetration depth with the XY 
exponent, in contradiction to the results obtained from the Ginzburg-Landau theory for the 
superconducting order-parameter. Possible reasons for this difference are discussed. 
PACS: 74.40.+k, 64.60.Ak 



Duality transformations in statistical mechanics have led to important insights into 
the nature of phase transitions in various systems, ever since the discovery of self-duality 
in two-dimensional (2D) Ising model by Kramers and Wannier This concept has also 
proved important in studies of the phase transition in 2D XY model, where dual partition 
function factorizes into spin- wave and vortex contributions 0, and therefore provides a 
natural framework for description of the Kosterlitz-Thouless transition. Similarly, the phase 
transition in 3D XY model can be understood within a dual approach, where the basic 
objects, instead of vortices, now become vortex loops. Critical temperature is identified as 
the point of proliferation of vortex loops through the system. The high-temperature phase 
of the original XY model corresponds to the low-temperature phase of the dual model, in 
which the disorder field that describes vortex loops has condensed 0|. Critical exponent 
for the correlation length calculated within the dual approach |^ is in very good agreement 
with the standard value, confirming the correctness of the physical picture of condensation 
of vortex loops at the transition in 3D XY model. 

Critical behavior of superconductors is described by the Ginzburg-Landau theory for a 
complex scalar field which. Cooper pairs being charged, is coupled to a fluctuating gauge po- 
tential. Nature of the phase transition in this system has long been a matter of debate: while 
close to four dimensions the transition is always first-order |P, P], in three dimensions, it is 
believed that it can be both first and second order, depending on the value of the Ginzburg- 
Landau parameter n 0, 0, [|10]. While some consensus has emerged on the overall 
topology of the phase diagram in the 3D case, the precise nature of the continuous transition 
in the theory is still controversial. Recent renormalization group analysis of the Ginzburg- 
Landau theory for superconducting order parameter yields a novel critical behavior, distinct 
from the one of the neutral superfiuid 0, |Ty]. Monte Carlo and analytical studies of the 



lattice version of the theory suggest, on the other hand, that the transition still belongs to 
the XY universality class ||1T|, [j^. As a step towards a better understanding of these issues, 
in this paper the theory for 3D lattice superconductor is revisited. Under a certain plausi- 
ble, but unproven assumption, the dual partition function of lattice superconductor may be 
expressed in a form that has a simple continuum limit. Based on different arguments, the 
same form of the continuum dual theory was previously proposed in ref. 12, and studied 
in ref. 13. However, there are some important differences between the assumptions and 
the conclusions of the present work and those of the previous study |T^. The continuum 
dual model is analyzed by renormalization group in 3D and the transition described by this 
theory is confirmed to be in the universality class of 3D XY model. The stiffness constant 
proportional to the superfiuid density vanishes as a power-law close to the transition, with a 
universal amplitude given by the fixed point value of the dual charge. In combination with 
the result on the anomalous dimension of the gauge field at the stable fixed point in the 
original theory [Q, the continuum dual theory suggests that the magnetic field penetration 



depth diverges with the power v^y — 2/3. This agrees with the result of Peskin |T^, but 
not with the conclusions of ref. 13. It also disagrees with the exponent derived from the 
Ginzburg-Landau theory 0, |T^. Possible reasons for this difference are considered. 

The partition function for lattice superconductor (3D lattice electrodynamics) is defined 
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by H, 

Z = |^'^exp(iroEcos(ft"+A-^n-A„,)-^5:(Vx (1) 

where the sums run over sites and hnks of a 3D quadratic lattice, fi = x,y,z, nfi denotes 
an oriented hnk between n and n + fi sites, integrals over gauge potential in the transverse 
gauge V ■ A = run over whole real axis, integrals over phases are from — tt to tt, and V x A 
and V ■ A are the lattice curl and the lattice divergence, respectively, e is the charge of a 
Cooper pair and Kq is the stiffness constant proportional to the superfluid density in the 
Meissner phase. The lattice superconductor corresponds to the extreme type-II limit of the 
Ginzburg-Landau model where the phase transition between normal and superconducting 
phases is expected to be continuous , [|10[ . In the Villain approximation 0] the phases can 
be integrated out exactly, and the theory may be expressed in terms of integer-valued link 
variables {m„^} coupled to the gauge field: 



where primes on the integral and on the sum, as throughout the paper, denote the conditions 
V ■ v4 = V ■ m = 0. The second condition forces the integer variables {mn^} to form 
closed loops. For e = the gauge field decouples, and one is left with the standard loop 
representation of the XY model, with the continuous phase transition at some Kq = K^. 
If Kq = oo, the sum over integers {mn^} can be performed, and the gauge-field becomes 
constrained to integer values. Remarkably, in this limit the lattice superconductor again 
maps onto XY model |M. Building on this observation, it has been argued that even for 



finite Kq the lattice superconductor should exhibit the "inverted" XY transition |Tl|. Here 
we note that the partition function can be rewritten as 

z = jJlT.^M-i^^T.('^^h)i^+iY.K,{wxAu-^^^^^ 

(3) 

so that the gauge field A can be easily integrated out. Since the field h is purely transverse, 
integer variables can be taken to be transverse as well, and 

[' ' 1 t 

Z = limt-,Q E exp(z27r E mn^,hn^, - ^ E(V ^ - E ^n/. - 2 E "^n/.)' (4) 

where we added a chemical potential term. Modulo Villain approximation, this is the same 



as 

f' f 1 1 
limt-,0 exp(- E cos(0„^ - 0^ - 2tiK^) - ^ E(^ ^ ~ 2" ^ 



When e = 0, the dual partition function (5) is identical to the Kq 00 limit of theory 
in (1). In this limit therefore, the transition is at = = 47r^frc, as originally found 
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by Peskin [Q. Apart from the limit t —>■ 0, the last expression is completely analogous to 
the partition function for the lattice superconductor, eq. 1, except that the dual gauge- field 
h is massive. Finite t adds a small chemical potential for the loop variables in the 

eq. 4. This suggests that the limit t — >■ is regular, in the sense that the nature of the 
phase transition is not changed by leaving t finite. This is the same assumption as the one 
made by the authors of ref. 11 in their argument for the inverted XY transition in lattice 
superconductor. Unfortunately, we have been unable to check directly the validity of this 
assumption. Rather, we will adopt it for a moment, and return to this question after we 
derive its consequences. For t ^ 0, the dual model of the lattice superconductor should be in 
the same universality class as the continuum theory for a complex scalar field representing 
vortex loops (the disorder field) minimally coupled to the massive dual gauge-field h (with 
V-K = 0) [0, H: 

L = [ rf^KV - /I)^(f)|2 + ^2|^(f)|2 + ^|^(f)|^ + IP + (V X hf]. (6) 

IJ? ~ (Tco — T) , Tco is the mean-field transition temperature, the quartic term describes the 
short-range repulsion between the vortex loops, and = c^Kq. 

Having a continuum version of the dual theory, the critical behavior can be studied by 
standard renormalization group methods (RG). We adopt the multiplicative field-theoretic 
renormalization and calculate the RG factors perturbatively and directly in 3D. The pertur- 
bation theory is used only to make the reasoning more explicit, and is not essential for the 
main conclusions. The renormalized dual theory is: 



/ Sr[Z^\{V - i-hm^ + m'2|vl/|2 + ^|v^|4 + \}? + ^(V X hf]. (7) 



Notice that the form of minimal coupling between and h is preserved and the same Ward 
identities hold as if h was massles. This further implies that the polarization of the field h is 
transverse and /i^ term does not get renormalized. To the leading order in two dimensionless 
coupling constants = bo/m and = (27r/e)^(/i|/m) we obtain: 



1+ 



247r' 
2go 



37r(l + (fih/m))' 



(9) 



6-^6o-AM_^(W^)\ (10) 
Stt m 2tt fih 

With the renormalized coupling constants defined as ^ = (27r/e)^(m^/m), = fi1/Zh, 
m = m' /Z]/'^ and b = b' /Z^, the one-loop /3-functions are: 

^^^lV^ = -9 + 7^9', (11) 
a log m 24:71 
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n - db h,^p 2 bg{2 + mh/m) , 1 .2 ..^n 

First, there are two unstable fixed points with vanishing dual charge Qc = 0: be = and 
be = Svr/S. Note that the last two terms in f3h depend not only on b and g, but on the third 
dimensionless combination of the coupling constants m/nih as well. The phase transition 
in the theory is achieved by tuning the temperature, that is by letting m — >■ 0, with the 
scaling of rrih determined by [3g. Close to the attractive fixed point with a finite dual charge 

{9 ~ 9c) = const ■ m^'a'^^''\ (13) 



so that 



^'^'^^ -ml = const ■ m^^^'a^^"^ + cjcm. (14) 
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Since (3g{gc) = 1 > 0, when m — the renormalized stiffness constant K = m|/e^ close to 
the transition approaches zero as 

Recalling that m = where is the correlation length of the disorder field, we see that 



the result is in accordance with the Josephson relation in 3D [|T6|. To the lowest order in 
dual charge, the universal amplitude is 

K 6 

hmm^o— = -. (16) 
m TT 

Now it becomes obvious that the last two terms in f3h vanish as the critical temperature 
is approached because mh/m — > 00, in spite of the scaling of dual charge towards a finite 
fixed point value. Consequently, the attractive fixed point with 7^ is located at the 
same be = Svr/S as the pure XY fixed point. As may have been expected, the fluctuations of 
the massive dual gauge-field do not influence the critical behavior of the theory, which is in 
the universality class of 3D XY model [|T2[, [|TH|. This is in sharp contrast with the effects 



of the massles gauge-potential in the Ginzburg-Landau theory [||, 0, ||T0[. We proved 



this statement here only to the lowest order, but clearly it must persist to all orders in 
perturbation theory. One could integrate out the dual gauge field h from the beginning 
and the only effect of that (apart from introducing some irrelevant terms) would be the 
modification of the bare quartic term for the disorder field. 

From the partition function for the lattice superconductor, eq. 1, it follows that the 
magnetic field penetration depth A scales with the stiffness constant as : 

where rjA is the anomalous dimension acquired by the gauge-field at the stable fixed point in 
the original theory. At the neutral, unstable fixed point with e = 0, t^a = 0. At the charged, 
attractive fixed point in the theory however, = 1 in 3D This may also be understood 



4 



by recalling that in the continuum limit the dimension of the charge [e^] = L~^, so that close 
to the stable fixed point scaling of the (V x A)^ term with length is changed by one power. 
From the relations (15) and (17) it follows that 

Aoce, (18) 

and the exponents for the correlation length of the disorder field and for the magnetic field 
penetration depth are the same. This result is independent of the perturbation theory. It 
was only necessary that (3g{gc) > 0, which is always true at a non-trivial zero of f3g. The 
continuum dual theory (6) thus implies that the exponent for the magnetic field penetration 
depth has the XY value u^y — 0.67. 

The same form of the continuum dual theory (6) has previously been considered in ref. 
13. However, by making some phenomenological assumptions for the coupling constants in 
the theory, these authors obtained an incorrect value for the exponent for the penetration 
depth. The essential difference with the present work is that the bare mass of the dual gauge 
field h was assumed to independently vanish close to the critical point (Tc — T) OC m^, 

X = l/i'xy The lowest order /3-function for the dual charge is then modified to read: 

Pg = (x-m-^f), (19) 

and the neutral {qc = 0) fixed point of the dual theory becomes infrared stable for any x > 1. 
Since Pg{0) = x — 1, eq. 13 then implies that 

~ m^. (20) 

The power x does not get changed by the fiuctuations of the disorder field if its assumed 
value is larger than unity. The authors of ref. 13 interpreted this observation as that the 
exponent for the penetration depth has the mean-field value 1/2. It is obvious however, that 
any other assumption for the power x would, as long as x > 1, lead to a correspondingly 
different exponent for the penetration depth. Furthermore, the interpretation of the field 
theory (6) as a continuum limit of the dual lattice superconductor makes it clear that nothing 
should be assumed for the mass fih- Close to the transition, the renormalized mass rrih is 
then driven to zero by fiuctuations of the disorder field, as a power law that follows from the 
theory. 

The behavior of the penetration depth derived here on the basis of the continuum dual 
theory is in disagreement with of results obtained from the Ginzburg-Landau theory for the 
superconducting order-parameter. Perturbative |^ and non-perturbative [|l^ renormaliza- 
tion group calculations in fixed dimension indicate that the penetration depth in Ginzburg- 
Landau model diverges with the exponent 0.50 < u < 0.62. The main point is that the 
exponent obtained from Ginzburg-Landau theory is different (and most likely smaller) than 
u^y. In fact, this is what one would naively expect, since the neutral, XY fixed point in 
the Ginzburg-Landau theory is unstable with respect to the charge, and small perturbations 
in that direction are attracted by the stable, charged fixed point with different exponents. 
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One conceivable reason for this disagreement is the approximate nature of the calculations 
performed in refs. 9 and 10. One should note however, that the non-perturbative method 
employed in ref. 10 leads to very accurate values for the exponents in several other prob- 
lems, where there is more information available. A more interesting possibility is that the 
condition t ^ which we were forced to relax in order to arrive at the continuum theory, 
might represent a singular limit in the theory |l^. In the eq. 4, for t = 0, the summation 



over integers {m„^} would force the dual gauge field h to take strictly integer values. Finite 
t relaxes this constraint, and h becomes distributed around integers with a finite width of 
distribution proportional to t. When viewed this way, it becomes less clear that finite and 
zero t necessarily lead to the same physics. More detailed numerical studies should be able 
to shed some light on the problem. 

The divergence of the penetration depth close to the superconducting transition is ex- 
perimentally measurable, at least in principle. In practice, unfortunately, this is made pro- 
hibitively difficult by smallness of the critical region for the gauge-field fiuctuations, which 
is rather narrow (even in the high-temperature superconductors) due to small value of the 
effective charge (Ginzburg-Landau parameter k, » 1). The presently accessible critical re- 
gion corresponds to the vicinity of the unstable neutral fixed point in the Ginzburg-Landau 
theory, at which the penetration depth diverges with the exponent I'xy/'^ (because rjA = 
at this fixed point) |jl8[. We should note here however, that there are other physical sys- 



tems, where the critical region for the gauge- field fiuctuations is much wider [|T9l, and 
which offer more hope for the experimental resolution of the question of the exponents at 
the charged critical point. 

In summary, the message of this paper is twofold. First, we studied the phase transition 
in 3D lattice superconductor within a continuum version of the dual theory which yields 
the transition to be in the XY universality class. The superfluid density vanishes as inverse 
of the correlation length of the disorder field. In contrast to the previous claims, dual 
charge is a relevant coupling whose fixed point value determines the universal amplitude 
for the scaling of the stiffness constant. This amplitude is calculated to the lowest order 
in perturbation theory. Finally, we have shown that within the continuum dual theory 
the magnetic field penetration depth diverges with the exponent u^y — 2/3. Second, it is 
pointed out that the critical behavior obtained on the basis of the continuum dual theory is in 
disagreement with the recent results from Ginzburg-Landau theory for the superconducting 
order-parameter. Possible reasons for this disagreement are discussed and experimental 
consequences are mentioned. 

Many useful discussions with Professors Z. Tesanovic and I. Affleck are greatfuUy ac- 
knowledged. This research has been supported by Natural Sciences and Engineering Research 
Council of Canada and Izaak Walton Killam foundation. 
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